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ABSTRACT 


Two consecutive three-day arcs of unsmoothed VHF Doppler tracking data from 
the Goddard Range and Range Rate System (GRARR) were used to determine 
lunar orbits for AIMP-E (Explorer 35). This study yielded magnitudes for 
Doppler residuals and estimates for trajectory errors as a function of five lunar 
gravity models. To a certain extent this analysis permits evaluation of the 
GRARR System tracking a lunar orbiting satellite. An additional purpose of this 
report is to relate orbit prediction quality as a function of various lunar gravity 
models. The trajectory differences were best minimized by the Field 15 
potential model which was derived from an analysis of orbital motion of Apollo 8. 
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COMPARISON OF LUNAR GRAVITY FIELDS USING GSFC 
RANGE-RATE DATA FROM AIMP-E 

1. INTRODUCTION 

The orbit of AIMP-E (launched on July 19, 1967) is markedly different from that 
of the other lunar satellites in the Apollo and Lunar Orbiter series. The orbit 
has a minimum height above the lunar surface of 790 kilometers with an eccen- 
tricity of 0,5739. Since AIMP-E has such a high orbit, lunar mass concentration 
effects are negligible and only the major terms in the lunar gravitational poten- 
tial and the perturbations due to the sun and the earth effect the AIMP-E orbit. 

The STADAN network is still tracking AIMP-E in lunar orbit using the VHF 
GRARR system. Some studies have been made to evaluate the lunar satellite 
tracking capability of the GRARR VHF system. Since future lunar orbiting satel- 
lites, such as Radio Astronomy Explorer (RAE-B), will also use VHF range-rate 
data, the following analysis also has a direct application. 

There are two practical methods of estimating the lunar satellite tracking capa- 
bility of Stadan VHF GRARR range-rate data. One method is to determine the 
magnitude of the Doppler frequency residuals in a converged solution for the 
AIMP-E orbit. Another approach is to determine the accuracy with which one can 
project or predict an orbit based on VHF Doppler observations. In this report 
trajectory evaluations were performed for each of five lunar potential models. 

The magnitudes of the position differences obtained in the various trajectory 
intercomparisons give an indication of the prediction capability of each gravity 
model. 

2. AIMP-E ORBITAL SOLUTIONS 

To study the orbit of AIMP-E, two consecutive three-day data arcs of VHF two- 

way Doppler data were obtained. The data were obtained by tracking stations 
operating at an uplink frequency of 148.26 MHz and a downlink frequency of 
12/13 the uplink frequency. The summary of contents for each data arc is listed 
in Table 1. Included in this table are the stop and start times of each pass. For 
each of the five gravity models, the two arcs of Doppler data were used to deter- 
mine an orbit. The numerical integration was performed with an Adams-Cowell 
twelfth-order integrator having a two-minute time step. The starting times for 
the differential correction on both data arcs are given in Figure 1, 

The five gravity fields which were used in the study were the tri -axial, R2, LI 
and R6 fields. Field 15, which was derived by Felsentreger, Murphy, Ryan and 
Salter, was also tried (reference 1). The coefficients of the five lunar potential 
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tricity of 0.5739. Since AIMP-E has such a high orbit, lunar mass concentration 
effects are negligible and only the major terms in the lunar gravitational poten- 
tial and the perturbations due to the sun and the earth effect the AIMP-E orbit. 

The STADAN network is still tracking AIMP-E in lunar orbit using the VHF 
GRARR system. Some studies have been made to evaluate the lunar satellite 
tracking capability of the GRARR VHF system. Since future lunar orbiting satel- 
lites, such as Radio Astronomy Explorer (ILAE-B), will also use VHF range- rate 
data, the following analysis also has a direct application. 

There are two practical methods of estimating the lunar satellite tracking capa- 
bility of Stadan VHF GRARR range-rate data. One method is to determine the 
magnitude of the Doppler frequency residuals in a converged solution for the 
AIMP-E orbit. Another approach is to determine the accuracy with which one can 
project or predict an orbit based on VHF Doppler observations. In this report 
trajectory evaluations were performed for each of five lunar potential models. 

The magnitudes of the position differences obtained in the various trajectory 
intercomparisons give an indication of the prediction capability of each gravity 
model. 

2. AIMP-E ORBITAL SOLUTIONS 

To study the orbit of AIMP-E, two consecutive three-day data arcs of VHF two- 

way Doppler data were obtained. The data were obtained by tracking stations 
operating at an uplink frequency of 148.26 MHz and a dowalink frequency of 
12/13 the uplink frequency. The summary of contents for each data arc is listed 
in Table 1. Included in this table are the stop and start times of each pass. For 
each of the five gravity models, the two arcs of Doppler data were used to deter- 
mine an orbit. The numerical integration was performed with an Adams-Cowell 
twelfth-order integrator having a two- minute time step. The starting times for 
the differential correction on both data arcs are given in Figure 1. 

The five gravity fields which were used in the study were the tri -axial, R2, LI 
and R6 fields. Field 15, which was derived by Felsentreger, Murphy, Ryan and 
Salter, was also tried (reference 1). The coefficients of the five lunar potential 
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models are listed in Table 2. The program input constants are listed in Table 
3. The geocentric coordinates of the VHF transmitting stations are given in 
Table 4. 


3. DOPPLER RESIDUAL MAGNITUDES 

The first gravity field used with AIMP-E Doppler data was the LI field. 

After four iterations of the LUNGFISH (Lunar Gravity Field in Spherical 

Harmonics) computer program, a converaged solution for the orbit was obtained 
using the second data stretch (August 23 to 26, 1968). The Doppler frequency 
residuals are plotted on Figures 2-5. The pass for station 52 (i.e., Carnarvon) 
has a duration of four minutes. The plotted residuals are bounded by an envelope. 
The envelope results because the range- rate data is recorded at a sampling 
frequency of one measurement per second which is shorter than the spin period 
of 2.41 seconds (reference 3). The envelope period of about 48 seconds contains 
an integral number of 20 spin periods. 

One method of gauging the quality of STADAN VHF Doppler data would be to com- 
pare the sum of squares and the sigma of the Doppler frequency residuals for 
each converged run. These two quantities are tabulated for each gravity model 
using both data arcs in Table 5. The values of sigma run from .4 to .5 Hz. For 
STADAN VHF data, one Hz of two-way Doppler frequency is equivalent to 1.1136 
meters/sec. The residuals, therefore, correspond to 40 to 50 cm/sec. In this 
study it was found that these residual magnitudes were of little significance in 
evaluating VHF range rate data. If all modeling errors are eliminated (ionosphere 
correction, satellite spin modulation and spin bias) the VHF GRARR system can 
achieve a resolution of 6 cm/sec using a sampling rate of 1 second (reference 5). 
The high noise level of 40 to 50 cm /sec obtained in the LUNGFISH runs is the 
result of AIMP spin cyclic modulation (in contrast to spin bias) effects. Spin 
modulation produces range-rate residuals that contain periodicities related to the 
spin of the satellite and the data sampling rate. The periodicities are caused by 
variations in the phase pattern of the antenna. In addition antenna rotation intro- 
duces a fixed bias in range rate measurement, (reference 3). The observation data 
was unsmoothed and data smoothing was left to the LUNGFISH orbit determination 
program. This procedure does not smooth spin bias. Thus the residuals appear 
somewhat high but this should not detract from the orbit computation accuracy. 

The magnitude of the residuals can be greatly reduced by fitting polynomials to 
the unsmoothed data prior to orbit computation. However, such data compression 
techniques were not employed during this analysis in order to preclude intro- 
ducing any bias due to smoothing. 

4. ORBIT COMPARISONS 

Another method of evaluating STADAN VHF GRARR data is to determine how 
accurately orbits can be predicted using this data. Such a study also permits 
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comparison of the predicting capability for the various lunar gravity n todels. 

For each gravity field, a propagated orbit (August 20 to 26) and a solution orbit 
(August 23 to 26) were obtained by differential corrections using the first and 
second data arcs respectively. The converged solution for the starting vector 
for each run is given in selenographic Kepler elements in Table 6. 

It is obvious that the starting vectors are almost identical for the tri- axial, R2 
and LI fields. The starting vectors for the R6 and Field 15 potential models 
are slightly different from those predicted by the previous fields. Since the 
starting vectors are so similar, one can expect the deviation of the solution 
orbit from the predicted orbit to be fairly similar for the five gravity fields. 
These results were expected because AIMP-E is in a high lunar orbit and the 
fields contain similar low-order and low-degree coefficients. 

The deviation of the solution orbit from the propagated orbit was transformed 
from selenocentric coordinates to radial, along-track and cross-track position 
differences. These trajectory differences are plotted on Figures 6-11. The 
graphs of the trajectory differences using each field exhibit a variety of forms. 
Most of the graphs have curves with a period equal to the orbital period of 
AIMP-E. Some curves appear to be periodic in nature, while others exhibit 
combinations of secular and periodic effects. For example, the position dif- 
ferences for the LI field and Field 15 vary in the following ways: The LI field 
has radial position differences that range from- 1.1 to + 1.3 km per orbital 
period at August 22, 19 hours GMT, 1968. Five orbits later, the radial position 
differences range from a minimum of -3 km to a maximum of +3.3 km per orbital 
period. The along-track error differences range from 2.1 to 4.0 km per orbital 
period at August 22, 19 hours GMT while five orbits later the maximum and 
minimum are 9.3 and 3.8 km respectively. The cross-track position differences 
appear to be periodic with maxima and minima of +12 and -41 km respectively . 
(The position differences of the R2 field are almost identical to those of the LI 
field). In comparison, Field 15 had a radial position difference curve which ap- 
pears to be periodic with an amplitude ranging from -1.1 to 1.5 km. The along- 
track position differences vary slowly from .4 to 6.0 km per orbital period at 
August 22, 19 hours GMT, 1968 and six orbits later vary from -.9 to 4.9 km per 
period. The cross-track differences appear to be periodic ranging from -44 to 
+14 km. Field 15 is seen to be superior in that over the three-day time interval, 
the radial and along-track position differences for Field 15 are almost periodic 
maintaining about the same maximum amplitude and form. For the other fields , 
the maximum amplitude of the radial and the along-track trajectory differences 
increase with each orbital revolution of AIMP-E. The largest position difference 
component for each of the five-fields is the cross-track component probably 
because range-rate supplies little information on the orientation of the orbital 
plane. 
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Field 15 does minimize the differences better than the other fields. However 
another factor in addition to gravity models seems to be very important in this 
study of trajectory differences. Since AIMP-E has a highly eccentric orbit, it 
spends a great deal of time near apolune. This causes certain problems in 
obtaining uniformly distributed data over the orbit. Due to satellite tracking 
schedules followed by the stations very few passes of data for AIMP-E last longer 
than five minutes. Since the orbital period is 11.52 hours, these passes form 

very short arcs. In the data processed for this report apolune was located at seleno- 
graphic longitude 80°E on August 19, 19 hours and at longitude 45° Eon August 22, 
16 hours. All of the data in each arc was taken when the satellite was 55 degrees 
on either side of apolune. If data were taken over the entire visible portion on the 
AIMP-E orbit rather than over just one quadrant, then the magnitude of the tra- 
jectory differences would be considerably reduced. This would require knowledge 
of the apolune and periiune times of AIMP-E at the tracking stations. 


5. SUMMARY 

This work accomplished two objectives. It served as an evaluation of STADAN 
VHF GRARR tracking performance in absence of modeling spacecraft spin and 
.ionospheric range-rate bias. If these quantities were modeled, the magnitude 
of the residuals of the original unsmoothed data would be greatly reduced. It is 
expected that the noise in the range-rate residuals would then be reduced from 
40 to 50 cm/sec and would approach the resolution of 6cm/sec (if a 1 second 
sampling rate were used) which is expected of the VHF GRARR system. These 
reduced residuals would be a measure of the accuracy that can be achieved. The 
second objective of the study was to indicate to what degree various lunar gravity 
fields influence orbit computation based on unsmoothed VHF GRARR data. The 
results of this study show that Field 15 (reference 1) minimizes the magnitude 
of the trajectory differences for AIMP-E better than the other four lunar gravity 
models using the August 20 to 26 (1968) AIMP-E range-rate data stretch. More 
uniform distribution of data over the orbit should further help to minimize the 
trajectory differences. 
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8. ABBREVIATIONS 

AIMP-E Anchored Interplanetary Monitoring Platform- E (Explorer 35) 

d Day 

GMT Greenwich Mean Time 

GRARR Goddard Range and Range Rate System 

hr or h Hour 

Hz Hertz, a cycle per second, unit of frequency 

km Kilometer 


LUNGFISH A computer program that was written by Computer Usage Co. 

for Langely S.F.C. to determine the gravity field of the moon 
from Lunar Orbiter range and range- rate data. At the pro- 
gram’s inception it handled only Deep Space Network data. At 
present it can also process MSFN and STADAN data. The mne- 
monics stand for Lunar Gravity Field in Spherical Harmonics. 




min or m 


Minute 


MHz 

R c 

RAE 

sec or s 

STADAN 

Me ' M © 1 M© 

VHF 


Megahertz 
Radius of the moon 

Radio Astronomy Explorer (Lunar Mission) 

Second 

Space Tracking and Data Acquisition Network 
Gravitational constants of the moon, earth and sun 
Very High Frequency — Frequency range of 30 to 300 MHz 


6 


Table 1 


Summary of AIMP-E Tracking Data (Two Arcs) 
DATA ARC 1 


Pass No. 

No. Obs. 

Start-Time 

End-Time 

Station 

Day 

H/M/S 

Day 

H/M/S 

1 

274 

233 

— 

024802 

233 

025301 

28 

2 

218 

233 

142005 

233 

142451 

27 

3 

291 

233 

160802 

233 

161304 

26 

4 

375 

233 

180621 

233 

181119 

28 

5 

292 

234 

013327 

234 

013822 

52 

6 

236 

234 

033219 

234 

033722 

52 

7 

285 

234 

052211 

234 

052657 

52 

8 

238 

234 

150338 

234 

150832 

26 

9 

295 

234 

170119 

234 

170617 

28 

10 

187 

235 

003103 

235 

1 

003618 

52 

11 

31 

235 

024415 

235 

024445 

52 

12 

319 

235 

120134 

235 

120704 

26 

13 

268 

235 

140111 

235 

140542 

27 

14 

284 

235 

160122 j 

235 

160617 

28 

15 

244 

235 

180948 

235 

181359 

27 

16 

300 

235 

233117 

235 

233616 
i 

28 


Data Arc 1 contains 16 passes and 4137 usable observations scattered over six 
11.52 hour orbits of AIMP-E 





Table 1 (Continued) 


DATA ARC 2 


Pass No. 

— i 

No. Obs, 

Start- Time 

End-Time 

Station 

Day 

H/M/S 

Day 

H/M/S 

17 

48 

236 

014747 

236 

014958 

52 

18 

308 

236 

033451 

236 

033958 

52 

19 

171 

236 

054155 

236 

054509 

52 

20 

271 

236 

130058 

236 

130544 

27 

21 

303 

236 

150514 

236 

151200 

26 

22 

184 

236 

171657 

236 

172000 

28 

23 

299 

236 

223136 

236 

223634 

28 

24 

375 

237 

003113 

237 

003612 

28 

25 

296 

237 

024239 

237 

024734 

52 

26 

287 

237 

043803 

237 

044249 

52 

27 

307 

237 

155113 

237 

155603 

27 

28 

268 

237 

212147 

237 

212642 

28 

29 

243 

237 

231755 

237 

232458 

28 

30 

300 

238 

012807 

238 

013306 

52 

31 

235 

238 

032027 

238 

032425 

52 

32 

276 

238 

145454 

238 

145953 

27 

33 

279 

238 

201629 

238 

202127 

27 

34 

287 

238 

221623 

238 

222121 

28 

35 

307 

239 

021807 

239 

022313 

52 


Data Arc 2 co ^ains 19 passes and 5044 usable observations scattered over six 
11.52 hour orbits of AIMP-E 
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Table 2 


/ 


y 

4 

5 

i 



Lunar Gravity Models 


Cor S 


71 

1 

Tri-Axial 

ii2 

LI 

H6 

Field 

N 

Field 

Field 

Field 

Field 

15 

C 

0 

w 

1 

* 

o j 

-207.1H07 

-207.108 

-207.108 

-210. 

-198.7651 

C 

:i 

0 1 

i 


21.0 

21.0 

16. 

- 75.66981 

C 

4 

0 






C 

2 

l 






C 

3 

l 


34. 

34. 

32. 

53.86133 

c 

4 

l 






c 

2 

2 

20.23974 

20.716 

20.716 

21. 

49.0906 

c 

3 

2 






c 

4 

2 






c 

3 

3 



2.583 



c 

4 

3 






c 

4 

4 





.1139595 

s 

2 

1 






s 

3 

1 




10. 

- 29.08456 

s 

4 

1 






s 

2 

2 





- 14.05312 

s 

3 

2 






s 

4 

2 






c 

kJ 

3 

3 






s 

4 

3 






s 

4 

4 





- .2008823 


* Multiple all coefficients by TO’ 6 
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Table 3 

List of Constants 



Table 4 

Geocentric Station Coordinates 


Transmitting 
Station NR 

R 

Radius 

(km.) 

0 

PHI 

Latitude 
1 11 

LAMBDA 
Longitude 
0 1 11 

22 

TAN 

0.6377320D 04 

-18 

54 

06.22 

47 

18 

12.56 

26 

ROS 

0.6371978D 04 

35 

00 

50.13 

277 

07 

26.89 

27 

SAN 

0.6372503D 04 

-32 

58 

30.90 

289 

19 

59.78 

28 

FRB 

0.6361008D 04 

64 

49 

25.96 

212 

29 

21.11 

52 

CRO 

0.6374450D 04 

-24 

45 

30.85 

113 

42 

55.06 












Sigma for Doppler Frequency Residuals 
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Selenographic Keplerian Elements at Epoch 





Figure 1. Greenwich mean times associated with each of the two data arcs and with the propagated and solution orbit. 
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Figure 1 1. Cross-track difference as a function of time for four gravity models. 
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